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Numerical simulation results are presented which suggest that a class of non-adiabatic 
rapid passage sweeps first realized experimentally in 1991 should be capable of im- 
plementing a set of quantum gates that is universal for one-qubit unitary operations 
and whose elements operate with error probabilities P e < 10 — 4 . The sweeps are non- 
composite and generate controllable quantum interference effects which allow the one- 
qubit gates produced to operate non-adiabatically while maintaining high accuracy. The 
simulations suggest that the one-qubit gates produced by these sweeps show promise as 
possible elements of a fault-tolerant scheme for quantum computing. 
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1 Introduction 

During the years 1997-1998 a number of researchers p]-[Z] showed that under appropriate 
circumstances a quantum computation of arbitrary duration could be carried out with ar- 
bitrarily small error probability in the presence of noise and imperfect quantum logic gates. 
The conditions that underlie this remarkable result are that: (1) computational data is pro- 
tected by a sufficiently layered concatenated quantum error correcting code; (2) fault-tolerant 
protocols for quantum computation are used; and (3) all quantum gates used in the compu- 
tation have error probabilities P e that fall below a value known as the accuracy threshold P a . 
One of the central challenges facing the field of quantum computing is determining how to 
implement quantum gates with error probabilities satisfying P e < P a . The accuracy thresh- 



a To whom correspondence should be sent. 

6 In this paper all gate error probabilities are per-operation. 



1 



old has been calculated for a number of simple noise models yielding results in the range 
10 -6 < P a < 10~ 3 . For many P a ~ 10 -4 has become a rough-and-ready working estimate 
for the threshold so that gates are anticipated to be approaching the accuracies needed for 
fault-tolerant quantum computing when P e < 10~ 4 . A number of universal sets of quantum 
gates have been found [5]-[I2] and so the problem of producing sufficiently accurate quantum 
gates has shifted to producing a sufficiently accurate universal set of such gates. One well- 
known universal set consists of the single-qubit Hadamard, phase, and 7r/8 gates together 
with the two-qubit controlled-NOT gate [13]. The single-qubit gates in this set are sufficient 
to construct any single-qubit unitary operation. 

In this paper numerical simulation results are presented which suggest that an existing 
class of non-adiabatic rapid passage sweeps [2] should be capable of implementing a set of 
quantum gates Si that is universal for one-qubit unitary operations. The one-qubit gates in 
Si are the Hadamard, NOT, V p , and Vl/g gates. The universality of Si for one-qubit unitary 
operations is established by noting that the Hadamard gate is an element, and demonstrating 
that the phase and ir/8 gates can be constructed from the V p , Vwg, and NOT gates. This is 
done in Section [3] For each of the gates in Si , sweep parameter values are presented which 
simulations indicate will yield gates that operate non-adiabatically and with error probabilities 
P e < 10~ 4 . This level of accuracy is a consequence of controllable quantum interference effects 
that are generated by these sweeps [15] . We explain the optimization procedure used to search 
for sweep parameter values that (when successful) yield this high degree of gate accuracy. 

The outline of this paper is as follows. In the following Section we summarize the neces- 
sary background associated with this class of non-adiabatic rapid passage sweeps; Section [3] 
presents our simulation results for the different gates; and Section |4] discusses these results, 
their relation to existing work in the literature, and also current challenges. 

2 Twisted Rapid Passage 

We consider a qubit that couples to an external control field F(t) through the Zeeman inter- 
action: 

H(t)=<r-F(t) , (1) 

where er are the Pauli matrices. The sweeps we will be interested in are a generalization of 
those used in adiabatic rapid passage (ARP). In ARP the field F(i) in the detector frame [16] 
is inverted over a time To such that F(<) = bx+atz. In an NMR realization of ARP, as seen in 
the lab frame, the detector frame rotates about the static magnetic field Bo z. In the detector 
frame, z is chosen to be parallel to the rotation axis. In the rotating wave approximation 
the rf-magnetic field B r / in the lab frame lies in the x-y plane and rotates about the static 
magnetic field. The detector frame is chosen to rotate with B,.j so that in this frame the rf 
field is static and its direction defines x: B r y — bx. The inversion time To is large compared 
to the inverse Larmor frequency Wq 1 (viz. adiabatic), though small compared to the thermal 
relaxation time Tth (viz. rapid). It provides a highly precise method for inverting the qubit 
Bloch vector s = (er), although the price paid for this precision is an adiabatic inversion rate. 
We are interested in a type of rapid passage in which the control field F(t) as seen in the 
detector frame is allowed to twist around in the x-y plane with azimuthal angle (f>(t) while 
simultaneously undergoing inversion along the z-axis: 

F(i) = 6cos0(i)x + 6sin0(t)y + atz . (2) 
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Here —To/2 < t < Tq/2 and y = z x x. Note that any pair of orthogonal unit vectors in the 
x—y plane can be used for x and y. Different choices simply alter the value of <p(t = 0). As 
will be explained shortly, interesting physical effects arise when the twist profile <f>(t) is chosen 
appropriately. This type of rapid passage is referred to as twisted rapid passage (TRP). 
The first experimental realization of TRP in 1991 by Zwanziger et. al. [13] carried out the 
inversion adiabatically with <p(t) = Bt 2 . Since then, non-adiabatic TRP has been studied with 
polynomial twist profile 4>(t) — (2/n)Bt n [15], and controllable quantum interference effects 
were found to arise for n > 3. Zwanziger et. al. [17 implemented non-adiabatic polynomial 
TRP with n = 3, 4 and observed the predicted interference effects. In the following subsection 
we briefly summarize how these quantum interferences arise and refer the reader to Ref. [15] 
for further discussion. 

2.1 Controllable Quantum Interference 

In the Zwanziger experiments [14] 117] , a TRP sweep is produced by sweeping the detector 
frequency linearly through resonance at the Larmor frequency ujq: (f>det(t) = ujq + (2at)/h. 
The frequency of the rf-field cf> r f is also swept through resonance in such a way that <f) r f (t) = 
<t>det{t) - <fi(t), where <j){t) = (2/n)Bt n is the TRP twist profile. Thus, 

<]) r f(t) = <t>det(t) - 4>{t) 

= c^ + — + . (3) 
At resonance <j) r f(t) = ^o- Inserting this condition into eq. ([3]), it follows that at resonance: 

at ~2i = °- (4) 

As shown in Ref. 15j, for polynomial twist <fi(t) = (2/n)Bt n with n > 3, eq. (U)) has n — 1 
roots, though only the real- valued roots correspond to resonance. The various possibilities are 
summarized in Table [1] We see that: (i) for B > a qubit always passes through resonance 

Table 1. Classification of regimes under which multiple qubit resonances occur for polynomial 
twist tj>(t) = (2/n)Bt n with n > 3. 



(a) n odd: 

(b) n even: 

(a) n odd: 

(b) n even: 



1. B > 

2 resonances at t = 

3 resonances at t = 

2. B < 



2 resonances at 
1 resonance at 



- and t ■ 
and t ■ 

and t : 




; (a/KB)»-2 
:±(a/RB)^ 

- (a/h\B\)^ 



multiple times during a single TRP sweep; (ii) for B < multiple resonances only occur 
when n is odd; and (iii) the time separating qubit resonances can be altered by variation 
of the sweep parameters B and a. Ref. [T5] showed that these multiple resonances have 
a strong influence on the qubit transition probability. It was shown that qubit transitions 
could be significantly enhanced or suppressed by small variation of the sweep parameters, and 
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hence of the time separating the resonances. Plots of the transition probability versus time 
suggested that the multiple resonances were producing quantum interference effects that could 
be controlled by variation of the TRP sweep parameters. In Ref. [TS] the qubit transition 
amplitude was calculated to all orders in the non-adiabatic coupling. The result found there 
can be re-expressed as the following diagrammatic series: 

Lower (upper) lines correspond to propagation in the negative (positive) energy level and the 
vertical lines correspond to transitions between the two energy levels. The calculation sums 
the probability amplitudes for all interfering alternatives |19j that allow the qubit to end up 
in the positive energy level at time t given that it was initially in the negative energy level. As 
we have seen, varying the TRP sweep parameters varies the time separating the resonances. 
This in turn changes the value of each diagram in eq. ([5]), and thus alters the interference 
between alternatives in this quantum superposition. Similar diagrammatic series can be 
worked out for the remaining 3 combinations of final and intial states. It is the sensitivity of 
the individual alternatives/diagrams to the time separation of the resonances that allows TRP 
to manipulate this quantum interference. Zwanziger et. al. |17j observed these interference 
effects in the transition probability using liquid state NMR and found quantitative agreement 
between theory and experiment. It is the link between the TRP sweep parameters and this 
quantum interference that we believe makes it possible for TRP to drive highly accurate 
single-qubit gates that operate non-adiabatically. The results presented in Section [3] for the 
different single-qubit gates are found by numerical simulation of the one-qubit Schrodinger 
equation. We next briefly describe how these simulations are done [15] . 

2.2 Simulation Protocol 

As is well-known, the Schrodinger dynamics implements a unitary transformation U(t,to) of 
the initial quantum state \tp(to)): 

m)) = u(t,t )m )) . (6) 

An n-qubit quantum gate implements a fixed unitary transformation U on n qubits. The 
unitary transformations Uh, Up, and Umot carried out by the one-qubit Hadamard, 

phase, 7r/8, and NOT gates are, respectively, 

Urr/8 = ( q e 2/4 J U NO T = ( 1 Q ^ ' ^ 

All matrices are in the representation spanned by the computational basis states |0) and |1) 
which are chosen to be eigenstates of a z : 

a z \i) = (-iy\i) (i = 0,l) . 

To determine the dynamical impact of TRP, the 1-qubit Schrodinger equation is sim- 
ulated numerically in the non-rotating frame in which the Hamiltonian H(t) is given by 
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eqs. (HJ) and It is found that the numerical stability of the simulation is enhanced 

if we expand the state \ip(t)) in the instantaneous energy eigenstates \E±(i)) for which 
H(t)\E±(t)) — E±(t)\E±(t)) . Because of the direct connection between these states and H(t), 
they carry substantial dynamical information, and a substantial portion of the dynamics due 
to H(t) can be accounted for by choosing this basis. This makes the task of determining the 
remaining dynamics using the Schrodinger equation much simpler and the simulation more 
stable. We thus write: 



!</;(*)> =S(t) exp 



do [e_ - 



T /2 

lit) exp 



\E-{t)) 



To/2 



d9 (E + - 



\E+(t)) 



(9) 



Here j±(t) are the adiabatic geometric phases [20] associated with the energy levels E±(t), 
respectively, and 

j±(t)=i{E ± (t)\f t \E ± (t)) . 

Substituting eq. (O into the Schrodinger equation leads to the equations of motion for S(t) 
and lit): 



where 



dS_ 
~dt 

dJ 
dt 



Sit) 
T(t) 



-T*(t) exp 
r(t) exp 



/ | d65(6) 

-To/2 



! i d6S(0) 

To/2 



lit) 

Sit) , 



(10) 



(E + it)\-\E_it)) 



[7+(*)-7-(*)] 



and r*(t) = — (E-it)\d/dt\E + it)) . The qubit is initially placed in one of the initial instan- 
taneous energy eigenstates |'0(— 7o/2)) = \E±(—Tq/2)) which fixes the initial condition for 
S(t) and lit) through eq. Q. It proves useful to recast eqs. (fl0|) in dimensionless form. To 
that end one introduces the dimensionless time r = (a/b)t, the dimensionless inversion rate 
A = h\a\/b 2 , and the dimensionless twist strength rj n = i%B / 'a)(b/ 'a) n ~ 2 . The connection 
between these dimensionless simulation parameters and the experimental sweep parameters 
is given in Section [3l It is straightforward to show that the resonances in Table Q] occur at 
US]: 

r = , (11) 

and 

^ r 1 

t = (sgnry n )<"- : 



(12) 



though only the real-valued solutions of eqs. pT|) and (fl^)) correspond to qubit resonances. 
The dimensionless version of eqs. (|10l) are the equations that are numerically integrated. The 
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simulations allow us to determine the actual unitary transformation U a produced by a specific 
assignment of the TRP sweep parameters To, a, b, B, and n. Section l2~4l will explain how the 
sweep parameters are iteratively modified so as to make U a approach a target gate U t as closely 
as possible. The iterative procedure searches for a sweep parameter set which minimizes (an 
upper bound for) the error probability P e for U a relative to U t - We next explain how P e and 
its upper bound are determined. 

2.3 Gate Error Probability 

The following argument is for an iV-dimcnsional Hilbert space, though N = 2 will be the 
case of interest in this paper. As in Section 12.21 let U a denote the actual unitary operation 
produced by a given set of TRP sweep parameters and Ut a target unitary operation we would 
like TRP to approximate as closely as possible. Introducing the operators D — U a — Ut and 
P — D^D, and the normalized state \if>), we define \tp a ) — U a \4>) and \4> t ) — U t \ip). Now 
choose an orthonormal basis \i) (i = 1, . . . , N) such that |1) = \tpt) and define the state 
via 

m = iv*> + i&> (13) 

= |l) + |^> ■ (14) 
Inserting \^) = ^,=i e d*) m t° eq. CU) g ives 

|^a> = (l + ei)|l)+£>|i) . (15) 

Since \ipt) = |1) is the target state, it is clear from eq. ([15]) that the error probability P e (tp) 
for U a (i. e. TRP) is 

P e W = ]T|e 4 | 2 . (16) 

We define the error probability P e for the TRP gate to be 

P e = max P e (?/>) . (17) 

m 

From eq. |[T5]). 

I&> = D\tl>) 

and 

= {i>\D^D\i>) 

= Trp^P , (18) 

where = \ip)(ip\. On the other hand, 



N 
i=l 

|2 



i=l 

= | ei | 2 +P e (V) • (19) 
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Combining eqs. (fT5)) and (TlT)|) gives 



PeW = (^l^)-leil 2 

< (6/>l&/>> = Tr P4>P ■ 

Since P = D^D is Hcrmitian it can be diagonalized: P = O^dO and d — diag(di, . . . , djy). 
Thus 

where = Op^,0^ L e ^ ^ = m ax((ii, . . . , cZat), then direct evaluation of the trace gives 

N 
AT 

< (p^) .. = Trp v , = d* , 

i=i 

where we have used that Tr'p^ = 1. Thus P e (ip) < d* for a/Z states l^). From eq. (fT7|) . it 
follows that 

P e <d» , (20) 

so that the largest eigenvalue d* of P is an upper bound for the gate error probability P e . 
Finally, notice that P = D^D is a positive operator so that di > for % = 1, . . . , N. Thus 
eZ* < Tr P and so 

P e < d* <Tr P . (21) 

Although Tr P need not be as tight an upper bound on P e as d* , it is much easier to calculate 
and so is more convenient than d* for use in the sweep optimization procedure to be described 
next. 

2-4 Sweep Optimization Procedure 

To find TRP sweep parameters that yield highly accurate non-adiabatic one-qubit gates we 
used the multi-dimensional downhill simplex method (21] to search for sweep parameters that 
minimize the upper bound Tr P for the gate error probability P e . Although we simulated 
a number of different types of polynomial twist, all data presented in Section [3] will be for 
quartic twist, </>4(r) = (i]4/2X)t 4 , which yielded the best results. The sweep parameters 
for quartic twist are (A, 774) which can be thought of as specifying a point in a 2-dimensional 
parameter space. For quartic twist, the downhill simplex method takes as input 3 sets of sweep 
parameters which specify the vertices of a simplex in the 2-dimensional parameter space. The 
dynamical effects of the TRP sweep associated with each vertex is found by numerically 
integrating the one-qubit Schrodinger equation as described in Section 12.21 The output of 
the integration is the unitary operation U a that a particular sweep applies. The desire is to 
iteratively improve U a so that it approximates as closely as possible a target unitary operation 
Ut- For each U a we determine P = (U a - Ut) f {U a - U t ) and evaluate TrP. The downhill 
simplex method then iteratively alters the simplex (i. e. one or more of its vertices) until 
sweep parameters are found that yield a local minimum of Tr P. Because this minimum is 
not global, some starting simplexes will give deeper minimums than others. Though there was 
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no gaurantee, it was hoped that a starting simplex could be found that yielded Tr P < 10~ 4 . 
Some trial and error in specifying the starting simplex was thus required, though for one-qubit 
gates, the trial and error procedure eventually proved successful and we present our results 
in the following Section. 

3 Simulation Results 

All results presented below are for quartic twist 



where r, A, and 774 are the dimensionless versions of time t, inversion rate a, and twist strength 
B (Section [272J) . For convenience, we re- write their definitions here: 



(- 



, h\a\ (hb 2 \ 



The parameter b was introduced in eq. ([2]) and is the rf field amplitude in an NMR realization 
of TRP [TTl [15]. All simulations were done with A > 1 corresponding to non-adiabatic 
inversion pHHI], and with r = aT /b = 80.000. 

The translation key connecting our dimensionless simulation parameters and the exper- 
imental sweep parameters used in the Zwanziger experiments [14[ I17j was given in the Ap- 
pendix of Ref. [15]. We re- write the formulas for quartic twist here for convenience. Note that 
Zwanziger's symbol B is here replaced by B to avoid confusion with our use of the symbol B 
in this paper to denote the twist strength. The translation formulas are: 

ui = ¥ (24) 



n 

aT 

h 



(25) 



B = =f- (26) 

" = < 27 > 

L0 1 ±Q 

* - Mi • (28 ' 

In the experiments of Ref. Q7]: uj x = 393Hz; T = 41.00ms; A = 50 000ffz; and B was 
calculated from eq. ([28]) with 774 varying over the range [4.50,4.70] x 10~ 4 . 
Note that Up and U^/s (see eqs. ([7]) and flSJ)) can be re-written as 

Up = e i7T/i U NO TVp (29) 

U n/S = e™/ 8 U N oTV„ /a , (30) 



where 



= ( 6 n ) > ( 32 ) 
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and Unot is given in eq. ([5]). As will be seen below, our simulations produced Vp and V^/g, 
from which Up and L^/g can be constructed using eqs. l[29|) and (|30|) . respectively. TRP is 
thus used to construct the set of gates Si — {Uh, Vp, V^/g, Unot} which is universal for one- 
qubit unitary gates. We stress that all gates in this set are produced using a non-composite 
TRP sweep (eq. ©). The different gates result from different choices for the TRP sweep 
parameters. For each one-qubit gate, we present our best-case results and show how gate 
performance is altered by small variations in the sweep parameters. 



Hadamard Gate 

The sweep parameters A = 5.8511 and 774 = 2.9280 x 10~ 4 produce the gate U a whose real 
and imaginary parts are: 



_ , 0.708581 0.705629 \ 

,! I 0.705629 -0.708581 J [ ' 

, _ , 0.380321 x 10^ 9 -0.144317 x 10~ 4 \ , , 

*m{U a ) - [ _ 144317 x 10 - 4 o.420313 xlO- 9 J' [ > 



For comparison, the real and imaginary parts of the target Hadamard gate Ut — Ur are: 



_ , 0.707107 0.707107 \ , . 

"" ( ;; ! ~ I 0.707107 -0.707107 J 

'»'('») = ( ) ' (36) 



From {/ a and Uh we find Tr P = 8.82 x 10 -6 and so the gate error probability satisfies 
P e < 8.82 x 10 -6 . Table [2] shows how gate performance varies when the sweep parameters 
are altered slightly. Of the two sweep parameters, 774 variation is seen to have the largest 



Table 2. Variation of Tr P for the Hadamard gate when the TRP sweep parameters are altered 
slightly from their best performance values. The columns to the left of center have r)4 = 2.9280 X 
10 — 4 and those to the right have A = 5.8511. 



m 


A 


TrP 




A 






TrP 




2.9280 x 10~ 4 


5.8510 


7.22 x 10" 


5 


5.8511 


2.9279 x 10" 


-4 


7.03 x 10 


-4 




5.8511 


8.82 x 10" 


6 




2.9280 x 10" 


-4 


8.82 x 10 


-6 




5.8512 


1.84 x 10" 


5 




2.9281 x 10" 


-4 


6.14 x 10 


-4 



impact on gate performance. This will turn out to be true for the other one-qubit gates as 
well. Although TRP can produce a Hadamard gate whose error probability falls below the 
accuracy threshold P a ~ 10~ 4 , it is clear from Table that the sweep parameters must be 
controlled to 5 significant figures to achieve this level of performance. See Section[5]for further 
discussion this point. 
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Vp Gate 

As noted above, the target gate here is Vp, and Up follows from eq. 1(291 
5.9750 and ?/ 4 = 3.8060 x 10~ 4 produce the gate U a : 

-0.627432 x 10~ 2 0.706181 
0.706181 0.627432 x 10~ 2 

-0.284521 x 10" 10 0.708004 
-0.708004 0.694222 x 10~ n 



eters A 

Re{U a ) = 

Im{U a ) = 
From eq. (|3"Tj). the target gate Vp is: 
Re(V P ) = 
Im(V P ) = 



0.707107 
0.707107 

0.707107 
-0.707107 



The sweep param- 

(37) 
(38) 

(39) 
(40) 



From U a and V P we find TrP = 8.20 x 10~ 5 and so P e < 8.20 x 10" 5 for this gate. Tabled 
shows how Tr P varies when 774 and A are varied slightly. Again gate performance is most 

Table 3. Variation of Tr P for the Vp gate when the TRP sweep parameters are altered slightly 
from their best performance values. The columns to the left of center have 7)4 = 3.8060 X 10~ 4 
and those to the right have A = 5.9750. 



m 


A 


TrP 




A 






TrP 




3.8060 x 10~ 4 


5.9749 


1.56 x 10' 


4 


5.9750 


3.8059 x 10' 


-4 


2.29 x lO - 


:! 




5.9750 


8.20 x 10- 


5 




3.8060 x 10" 


-4 


8.20 x 10" 


5 




5.9751 


1.43 x 10" 


4 




3.8061 x 10" 


-4 


1.88 x 10" 


3 



sensitive to variation of 774, and the sweep parameters must be controlled to 5 significant figures 
for performance to surpass the accuracy threshold. The latter point is discussed further in 
Section |H 

Vx/g Gate 

From eq. ([30]) . U n /g is found from V^/ 8 and Umot- The target gate this time is K-/8- F° r 
A = 6.0150 and r/ 4 = 8.1464 x 10~ 4 TRP produced the gate U a : 



Re{U a ) = 

Im(U a ) = 
From eq. (|32|). the target gate V^ 
ReiV^/s) 
ImiV^/s) 



0.101927 x 10- 2 0.925307 
0.925307 -0.101927 x 10- 2 

-0.960223 x 10~ 10 



-0.379218 





0.923880 



0.379218 
0.184961 x 10~ 10 



0.923880 






-0.382683 



0.382683 




(41) 
(42) 

(43) 
(44) 



These matrices give Tr P — 3.03 x 10 5 and so for this gate P e < 3.03 x 10 -5 . Table g] shows 
how gate performance varies when the sweep parameters are altered slightly. As with the 
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Table 4. Variation of Tr P for the V n i% gate when the TRP sweep parameters are altered slightly 
from their best performance values. The columns to the left of center have 774 = 8.1464 X 10~ 4 
and those to the right have A = 6.0150. 





A 


TrP 




A 


»?4 




TrP 




8.1464 x 1CT 4 


6.0149 


1.30 x 10- 


•i 


6.0150 


8.1463 x 10" 


-4 


1.77 x 10 


-3 




6.0150 


3.03 x 10- 


5 




8.1464 x 10" 


-4 


3.03 x 10 


-5 




6.0151 


2.18 x 1(T 


3 




8.1465 x 10" 


-4 


2.77 x 10 


-3 



previous two gates, performance is most sensitive to variation of 774, and the sweep parameters 
must be controllable to 5 significant figures (see Section rjj. 

NOT Gate 

Finally, we examine U NO t- For A = 7.3205 and 774 = 2.9277 x 10~ 4 TRP produced the gate 
U a : 



, M _ / 0.235039 x 10~ 2 
He[Ua) ~ \ 0.999997 

j- (tt \ _ ( -0.323648 x lO" 10 
~ ^ 0.115150 x 10- 4 

For comparison, Umot is (eq. dHJ): 

ReQJNor) = 
Im(U NO T) = 



0.999997 
0.235039 x 10~ 2 

-0.115151 x 10" 4 
0.271006 x 10~ 10 



1 

1 







(45) 
(46) 

(47) 
(48) 



These matrices yield TrP = 1.10 x 10~ 5 and so P e < 1.10 x 10~ 5 . Table [5] shows how TrP 
varies with small variation of the sweep parameters. As with the other gates, performance is 

Table 5. Variation of Tr P for the NOT gate when the TRP sweep parameters are altered slightly 
from their best performance values. The columns to the left of center have 774 = 2.9277 X 10 -4 
and those to the right have A = 7.3205. 



774 


A 


TrP 




A 


V4 




TrP 




2.9277 x 10~ 4 


7.3204 


1.12 x 10" 


5 


7.3205 


2.9276 x 10" 


-4 


1.23 x 10" 


3 




7.3205 


1.10 x 10" 


5 




2.9277 x 10" 


-4 


1.10 x 10" 


5 




7.3206 


1.22 x 10- 


5 




2.9278 x 10" 


-4 


1.23 x 10" 


3 



most sensitive to variation in 774, and sweep parameters must be controllable to 5 significant 
figures for the gate error probability P e to fall below the accuracy threshold P a ~ 10 -4 (see 
Section gj). 

4 Discussion 

In this paper we have presented numerical simulation results which suggest that TRP sweeps 
should be capable of producing a set of quantum gates that is universal for one-qubit unitary 
operations. We also showed that sweep parameters can be found which the simulations 
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indicate will yield gates that operate non-adiabatically and with error probabilities satisfying 
P e < 1CT 4 . To achieve this de gree of accuracy, however, the sweep parameters must be 
controllable to high precision (5 significant figures). This raises the question of whether such 
precision is possible with current technology. In the NMR realization of TRP [14] [17] the 
inversion time T was of order 10~ 2 s, while the spectrometer waveform resolution allowed 
the rf- and detector-phases to be specified in time steps of order 10~ 7 s. Thus To can be 
determined to 1 part in 10 5 . By using shimming and sample rotation the uncertainty in the 
Larmor frequency ujq (which is used in eqs. (|53j) and (|54j) ) can be reduced to 15 Hz, while 
its value is 500MHz. This corresponds to a relative error of Awo/^o ~ 10~ 7 . It is also 
possible to use tt or tt/2 pulses to calibrate the rf field strength lo\ down to a relative error of 
Awi/wi ~ 10~ 4 . Thus many of the TRP sweep parameters are already at or near the level of 
precision needed to make high-fidelity one-qubit gates. Still, it is clear that further theoretical 
work is needed to find ways to make the gate error probability a more slowly varying function 
of the TRP sweep parameters. Recall that TRP sweeps are non-composite. It is an interesting 
open question whether composite sweeps that interlay TRP with different types of pulses can 
lead to more robust gate performance. We intend to examine this question in our next set 
of simulations. Having discussed current challenges, it is worth stressing that these sweeps 
show genuine potential for producing high-fidelity non-adiabatic one-qubit gates. Further 
work to try to develop this potential seems warranted. Although other approaches exist for 
making one-qubit gates (e. g. Ref. [23]), in a field faced with as many technical challenges as 
quantum computing, it is advantageous to have multiple ways to accomplish important tasks. 
It is hoped that with further development TRP gates may provide an approach to making 
high-fidelity non-adiabatic quantum gates. TRP sweeps also provide a concrete example of 
how quantum effects can be used to enhance our control of a quantum system. Further study 
of these sweeps also seems worthwhile as a question of basic physics. 

Atomic Physics 

The following scenario is inspired by the NMR realization of TRP p~4] [17] . Consider electric 
dipole transitions between a pair of atomic energy eigenstates |+) and |— ) of the Hamiltonian 
H a with respective energies E± = ±eo/2. Transition between these two states is caused by 
an applied electric field E a (i) = 2E\ cos cf> a (t) x which couples to the atom's electric dipole 
moment d = er. In the lab frame, the two-level Hamiltonian Ti(t) in the rotating wave 
approximation is |22| : 



where Hloq = eo arL d = d x E\. Transformation to the detector frame [141116) is done using 
the unitary operator U(t) = exp[— (i/2)<fidet(t)<J z ] so that H — > H: 



H(t) 



a z + —— [ cos 4> a (t) <J X + sin cj> a (t) a y ] 




ata z + bcoscj) n (t)<7 x + b sin <p n (t)cr., 



(49) 
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where 



at = ~ ( (j) d et - wo ) (50) 
b = ^ (51) 

<£n(*) = 4>a~4>det , (52) 

and cj> n (t) = (2/n)Bt n is the twist profile for polynomial twist. Eq. (|49|) gives = er • F(i), 
where F(i) is the control field for TRP appearing in eq. ((2|). Integrating eq. ([50)1 gives 4>det(t) 
which can then be inserted into eq. ([52")) so that 

at 2 

<t>det(t) = -^+w i (53) 
n 

at 2 2 

4>a{t) = —+uj t + -Bt n . (54) 
n n 

We see that programming the generator that produces E a (i) so that the phase <fi a (t) is given 
by eq. (154)) causes a TRP sweep to be applied to the atom in the detector frame. Note that, 
to insure the two- level approximation is valid, the frequencies <p n {t) swept through by the 
TRP sweep should not include the resonance frequency of any other pair of atomic energy 
levels since this would drive unwanted dynamics not included in TL{t). 

Previous Work 

Recently, Morton et. al. [23] showed how to use composite pulses to produce high fidelity 
single-qubit operations in electron paramagnetic resonance. The composite pulses are based 
on the BB1 corrective sequence [24]. Along with observation of non-decay of Rabi oscilla- 
tions and suppression of secondary Fourier components in the spin echo decay envelope, they 
compared an improved Carr-Purcell pulse sequence (in which BB1 composite 7r-pulses replace 
ordinary 7r-pulses) with the Carr-Purcell-Meiboom-Gill sequence. From the decay of the echo 
produced by the improved Carr-Purcell sequence they inferred a fidelity for the BB1 7r-pulses 
of T = 0.9999. The authors noted that this fidelity is ultimately limited by pulse phase errors. 

The fidelity in Ref. [23] is T = (l/2)Re [Tr (U^U t )]. It is possible to relate our TrP 
upper bound on P e to this fidelity. Recalling that P = (U a — Ut) (U a — Ut), we have 

TrP - Tr(2-\ylU t + U}U a 
= 4-2Re[Tr(£/ a t [/ t )] 
= 4(1 -JO , 



and so 



T=\--TrP . (55) 
4 



Using the results from Section [3] for TrP in eq. ([55| . we can determine the fidelity for the 
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TRP gates: 



T H = 0.9999 98 

T Vp = 0.9999 80 

T v ^ /& = 0.9999 92 

Tnot = 0.9999 97 



(56) 
(57) 
(58) 
(59) 



Future Work 

(a) We are currently exploring whether TRP can be used to make a two-qubit gate that 

will complete the one-qubit gates considered here to give a set that: (i) is universal for 
quantum computation; and (ii) has all gates operating non-adiabatically with fidelities 
that yield P e < P a . A progress report on this work will be given elsewhere. 

(b) Development of an approximate analytical approach to TRP would be very useful. We are 

not aware of any general tractable analytical approach to non-adiabatic rapid passage 
that could be used to find good starting simplexes for the sweep optimization procedure. 
It is because of this that we followed the numerical approach described above. 

(c) Constructing a theory for the optimum twist profile <j)(t) for a given quantum gate would 

also be a valuable contribution. To date, quartic twist has worked best, though we do 
not presently have arguments explaining why it will produce better gates than the other 
examples of TRP that we have considered, or whether some other profile will work even 
better. 

(d) It would be especially interesting if the simulation results presented above could be tested 

experimentally. One possibility might be to use state tomography to measure the output 
density matrix p exp = U a \ipo){ipo \ U\ resulting from an initial state \4>o), for each of the 
TRP generated gates U a presented in Section [3j Associated with each sweep is a target 
gate Ut and a corresponding target density matrix p t = Ut^o) {"^oWl . Having measured 
p exp , evaluate the fidelity T(p exp ,p t ) [25]: 



Although this fidelity differs from the one considered in Ref. [23], one might naively 
anticipate that they are of comparable size. If so, then the experimentally determined 
fidelities should be close to the fidelities given in eqs. ([5^|) ~ ([5^]) . 

(e) The simulation results presented in this paper are for an isolated qubit interacting with a 
noiseless TRP sweep. Although this scenario might appear idealized, it seemed sensible 
to see what kind of performance was possible using TRP under the best possible con- 
ditions. One important extension would be to allow the TRP sweeps to include a noise 
component. To the extent that this noise leads to dephasing, TRP gate performance is 
expected to deteriorate once the dephasing time is of order the TRP inversion time To. 
Under these conditions, the qubit dynamics begins to loose its temporal phase coher- 
ence, and the quantum interference between alternatives begins to disappear. It would 
be worthwhile to consider simple noise models to study the sensitivity of TRP gate per- 
formance to parameters such as noise power and noise correlation time. A study along 




(60) 
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these lines was done for the quantum adiabatic search algorithm in Ref. [55]. Phase 
decoherence resulting from interaction of the target qubit with environmental qubits is 
another source of concern. As with the case of noise above, should the decoherence time 
be of order To, TRP gate performance is expected to suffer. Follow-up work that sheds 
light on how this performance cross-over occurs would be valuable. 
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